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The iron chalcogenides FeTe and FeSe belong to the family of iron-based superconductors. We study 
the magnetism in these compounds in the normal state using the ab initio downfolding scheme developed 
for strongly correlated electron systems. In deriving ab initio low-energy effective models, we employ the 
constrained GW method to eliminate the double counting of electron correlations originating from the 
exchange correlations already taken into account in the density functional theory. By solving the derived ab 
initio effective models, we reveal that the elimination of the double counting is important in reproducing the 
bicollinear antiferromagnetic order in FeTe, as is observed in experiments. We also show that the elimination 
of the double counting induces a unique degeneracy of several magnetic orders in FeSe, which may explain 
the absence of the magnetic ordering. We discuss the relationship between the degeneracy and the recently 
found puzzling phenomena in FeSe as well as the magnetic ordering found under pressure. 


After the discovery of iron-based superconductors, 
many families of the iron-based superconductors have 
been found. 1-3 ) Among them, iron chalcogenides such as 
FeTe and FeSe (11-type) are the simplest iron-based su¬ 
perconductors, which have no blocking insulating layer. 
Despite the simple lattice structure of the 11-type com¬ 
pounds, many characteristic behaviors have been ob¬ 
served in them, which have attracted much interest. For 
example, FeTe shows an exceptional bicollinear antifer¬ 
romagnetic (AFB) order with an ordered moment of 
~ 2.0-2.25/xb, 4 ’ 5 ^ although most of the iron-based su¬ 
perconductors show a stripe antiferromagnetic (AFS) or¬ 
der. Although the AFB order itself is reproduced by sev¬ 
eral calculations based on the density functional theory 
(DFT), 6-9 ) the microscopic origin of the AFB order and 
the roles of electron correlations are not fully clarified 
yet. 

In FeSe, although calculation based on the DFT pre¬ 
dicts that the AFS order is the ground state, 8 ) no clear 
signature of magnetic order is observed in experiments, 
and the superconducting phase appears below T c ~8 
K. 10 ) Note that a structural phase transition without 
magnetic order occurs at T s ^90 K. The absence of the 
magnetic order despite its large local moment 11 ) remains 
a puzzle in FeSe. Recent experimental result 12-14 ) sup¬ 
ports the notion that FeSe has extremely small Fermi 
surfaces and that the resultant superconductivity is lo¬ 
cated in the Bose-Einstein-condensate (BEC) region. 
The origin of these exotic phenomena also remains a 
puzzle in FeSe. Furthermore, it is reported that a sin¬ 
gle layer of FeSe on SrTiC >3 shows high-temperature su¬ 
perconductivity, whose T c is higher than 50 K. 15 ) Thus, 
clarifying these remarkable behaviors in iron chalco¬ 
genides is important for understanding the nature of 
high-temperature superconductivity in iron-based super¬ 
conductors. 

Many theoretical and experimental studies indicate 
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that the 11-type compounds are located in the strong¬ 
coupling region. 11-14,16-20 ) In particular, in recent ex¬ 
perimental studies of a single crystal of FeSe, 12-14 ) it has 
been shown that the experimental Fermi surfaces can¬ 
not be explained by the conventional DFT calculations. 
In addition, an ab initio evaluation of the interaction 
parameters suggests that the 11-type compounds have 
relatively stronger correlations than other iron-based su¬ 
perconductors such as BaFe 2 As 2 and LaFeAsO. 17 ) In the 
strong-coupling region, it is necessary to consider the 
effects of electron correlations beyond the conventional 
DFT. Even when the strong electron correlation is taken 
into account, the above recently studied Fermi surface 
does not appear consistent with the calculated Fermi sur¬ 
face. 18,21 ) This result is a puzzling aspect of FeSe. 

A hybrid method (ab initio downfolding scheme) 22 ) of 
the DFT and precise model calculation is a promising 
way of treating electron correlations beyond the stan¬ 
dard DFT: In this scheme, we first obtain the global 
band structure based on the DFT and evaluate the inter¬ 
action parameters in the low-energy effective model by 
constrained random-phase approximation (cRPA). 23,24 ) 
Then, we solve the low-energy effective model using pre¬ 
cise low-energy solvers such as the many-variable varia¬ 
tional Monte Carlo (mVMC) method. 25 ) Several appli¬ 
cations to real materials show the validity and accuracy 
of this scheme 26-28 ). However, in this scheme, there is 
a fundamental problem in that the electron correlations 
are doubly counted, i.e., the electron correlations are con¬ 
sidered in both the band calculations and the model cal¬ 
culations. Although this double counting is assumed to 
be small in the conventional downfolding scheme, in the 
strong-coupling regime where iron chalcogenides are ex¬ 
pected to be located, the elimination of double counting 
may become important to correctly reproduce electronic 
structures. 

In this Letter, by employing the state-of-the-art ab 
initio method, we clarify the magnetic properties in iron 
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Fig. 1. (Color online) Electronic band structures of the Fe 3d 
MLWFs of FeTe in the LDA [(red) solid line] and the cGW-SIC 
[(green) dashed line]. The Fermi energy is set to zero. 


chalcogenides. We properly eliminate the double count¬ 
ing of the electron correlations and derive the low-energy 
effective models in an ab initio way. 29 ) Then, we solve 
the low-energy effective models using the mVMC. As a 
result, we show that the elimination of the double count¬ 
ing plays a crucial role in stabilizing the AFB in FeTe. 
We also solve the low-energy effective models for FeSe in 
addition to FeTe, and find that several magnetic orders 
are energetically nearly degenerate. This degeneracy may 
result in the absence of the magnetic order in FeSe un¬ 
der competitions and frustrations emerging in the region 
close to the phase boundaries. 

Our low-energy effective models for iron-based super¬ 
conductors are given as follows: 


4 = EE - Rj)4nAm„ + ? EE 


ij mncrp 


Umncrp(Ri ^j^^ J ina^\rnp^ J 3 rn P^ n(J Jmnap^Ri Rj ) 


{4na4mpdi n pdjma- + 4na4npdjmpdjma) f j (1) 


tmna(R) = (<£m0 \R X \(j> nR ) , (2) 

where t x represents the transfer integral of the maxi¬ 
mally localized Wannier functions (MLWFs). 30,31 ) Here, 
X denotes the methods of deriving the one-body part of 
the Hamiltonian, i.e., X=LDA (cGW-SIC) means that 
the local density approximation (constrained GW-self- 
interaction correction) is used for evaluating transfer in¬ 
tegrals. We will show details of the cGW-SIC later. The 
corresponding band structures near the Fermi level are 
shown in Fig. 1. Here, d\ na ( di ncr ) is a creation (anni¬ 
hilation) operator of an electron with spin a in the nth 
MLWF having strong Fe 3 d characters (with n = 1-5 
being the orbital index) centered at R;. 17 ) In this work, 
we take the X- and Y-axes for the MLWFs index along 
the Fe-Te/Se directions. Note that the 12 x 12 x 6 fc-mesh 
is employed in the LDA calculation, and the 6x6x3 


fc-mesh is employed in the cRPA and cGW-SIC calcu¬ 
lations. The interaction parameters U and J denote the 
Coulomb interactions and exchange interactions, respec¬ 
tively. These interaction parameters are calculated by the 
cRPA 22,23 ) and dimensional downfolding. 27,32,33 ) For de¬ 
tails of the interaction parameters and hopping param¬ 
eters, see Ref. 34. The off-site interactions (i ^ j ) are 
omitted in this calculation because they are less than 
one-quarter of the on-site interaction parameter. Double 
counting in the Hartree term is subtracted as we have 
mentioned in our previous paper. 26 ) 

For the one-body part of the Hamiltonian, in addi¬ 
tion to the standard LDA calculations, the constrained 
GW (cGW) scheme was employed, 29,35 ) which can elim¬ 
inate the double counting of the self-energy originating 
from the exchange-correlation energy. In this scheme, we 
first subtract the exchange-correlation energy in the LDA 
calculations and replace it with the GW self-energy cal¬ 
culated by excluding the low-energy contribution, since 
the low-energy part of the self-energy will be considered 
later when the obtained model is solved. In this scheme, 
since we do not consider the low-energy part of the GW 
self-energy in deriving the low-energy effective model, the 
double counting of exchange correlations does not occur. 

The cGW self-energy does not include the self¬ 
interaction (SI) originating from the low-energy degree of 
freedom, which should be canceled out by the opposite- 
sign SI in the Hartree term. In the nondegenerate multi¬ 
band systems, the imbalance of the SI exerts crucial ef¬ 
fects on the occupation number of each orbital and the 
magnetism. We therefore introduce the self-interaction 
correction (SIC) originating from the low-energy degrees 
of freedom U on ^ SIte n^ DA /2 into the cGW scheme, where 
U on-site j g ^j ie on _ s ite effective interaction and n ^ BA /2 is 
the occupation number of the low-energy degrees of free¬ 
dom of the up or down spin in the LDA. The one-body 
part with the SIC is given by 


^cGW-SIC 


^cGW 


_ ^cGW ^yon-site ^1.1)A 


( 3 ) 


where R cGW is the static one-body Hamiltonian in the 
cGW. 29 ) The renormalization factor Z cGW is needed to 
renormalize the frequency-dependent part of the interac¬ 
tion of the effective model. 29 ) 

Figure 1 shows the band structures of FeTe in the 
LDA and cGW-SIC, which correspond to the one-body 
part of the low-energy effective model in Eq. (1). We 
note that the band structure of the cGW-SIC model for 
the noninteracting case is a virtual one, which cannot 
be directly observed in experiments. It is expected that 
correct band structures are reproduced after solving the 
low-energy effective models. The main difference between 
the LDA and cGW-SIC models is the on-site potential of 
the 3 Z 2 — R 2 orbital. In the cGW-SIC model, the on-site 
potential of the 3 Z 2 — R 2 orbital decreases and the occu¬ 
pation of the 3 Z 2 — R 2 orbital increases. As we will show 
later, the stability of the magnetic order largely depends 
on the on-site potential of the 3 Z 2 — R 2 orbital. 

To solve the low-energy effective models, we employ 
the mVMC, which can treat the electron correlations se- 
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Fig. 2. (Color online) AF patterns calculated in this work: (a) 
AFG, (b) AFB, (c) AFS, and (d) AFH. (e) Groundstate energy 
per site of AF in mVMC extrapolated to the thermodynamic limit 
for FeTe with t LDA and with t cGW ~ slG (in meV). The energy E 
is extrapolated to the thermodynamic limit Eoo by employing the 
scaling E = E 00 + aN s 3 / 228 > 36 ) with a constant a for the system 
size N s . 


riously. 25 ) Our variational wavefunction is defined as 

I i>) = T 3 V G C S =° Impair), (4) 

where Impair) is a generalized pairing wave function, and 
C s ~°, Vg, and V.j are the total spin quantum-number 
projection, Gutzwiller factors, and Jastrow factors, re¬ 
spectively. 25, 26 i For more computational details of iron- 
based superconductors, readers are referred to Ref. 26. 

In the mVMC, we examine the stability of four possible 
AF orders, i.e., G-type AF state (AFG), AFS, AFB, and 
half-collinear AF state (AFH), whose schematic pictures 
are shown in Figs. 2(a)-2(d), respectively. We mainly 
show the results of N s = 8 x 8 sites with periodic bound¬ 
ary conditions; the systems with 4x4 and 4x8 sites 
are also calculated for the extrapolation to the thermo¬ 
dynamic limit. 

By solving the low-energy effective models, we show 
how the elimination of double counting affects the sta¬ 
bility of magnetic orders in FeTe. As shown in Figs. 2(a)- 
2 (d), we examine four different magnetic orders, which 
are potential candidates of the order. In the LDA model, 
we find that the AFS is the most stable magnetically 
ordered phase [Fig. 2(e)]. In contrast, the AFB phase 
becomes the ground state in the cGW-SIC calculation 
that properly eliminates the double counting of the self¬ 
energy. This result indicates the importance of the elim¬ 
ination of the double counting in the strong-coupling 
regime. 

In Table I, we show the orbital resolved occupation 


Table I. Occupation numbers of the Wannier functions of the 
LDA and cGW-SIC models in the mVMC, where the sum of the 
occupancy is 6. Error bars are negligible (< 10 3 ). 


FeSe (AFH) 

XY 

YZ 

3Z 2 - R 2 

YX 

X 2 -Y 2 

LDA 

1.14 

1.36 

1.02 

1.32 

1.16 

cGW-SIC 

1.27 

1.04 

1.64 

1.04 

1.01 

FeTe (AFB) 

XY 

YZ 

3Z 2 - R 2 

YX 

X 2 — Y 2 

LDA 

1.40 

1.12 

1.12 

1.31 

1.05 

cGW-SIC 

1.44 

1.05 

1.42 

1.07 

1.02 


number for both the LDA model and the cGW-SIC 
model. Compared with that in the LDA model, the oc¬ 
cupation number of the YZ/ZX orbitals becomes close 
to half filling ( riyz/zx ~ 1) in addition to the X 2 — Y 2 
orbital. This is caused by the lowering of the on-site po¬ 
tential of the 3 Z 2 — R 2 orbital in the cGW-SIC model. 
The localization of the X 2 — Y 2 orbital is frequently re¬ 
ferred to as the orbital selective Mottness in the litera¬ 
ture. 18 . 26 > 2 b 37-39) i n addition, our result indicates that 
the YZ/ZX orbitals also show orbital-selective Mott in¬ 
sulating behavior in FeTe as well as in FeSe. The or¬ 
bital selective Mott insulating behaviors are evidenced by 
small double occupancies d m = (n m fn m i) = 0.1 for m = 
YZ/ZX and = 0.04 for m = X 2 — Y 2 and small orbital- 
resolved charge compressibilities Xcm = dn m /d\i < 0.001 
(meV )” 1 for m = YZ/ZX and X 2 - Y 2 in FeTe (AFB) 
(fx = dE/dN is the uniform chemical potential applied to 
all the orbitals, where E is the total energy and N is the 
total number of electrons). 4 °) Because the occupancies 
of the YZ/ZX orbitals are slightly different from half 
filling, the YZ/ZX orbitals can be regarded as doped 
orbital-selective Mott insulators in a strict sense, where 
the excess carriers from half filling make the system 
metallic . 41 ' 1 This additional localization may be the ori¬ 
gin of the characteristic behavior observed in iron chalco- 
genides. Note that a similar localization is observed in 
previous calculations for iron chalcogenides. 20,39 -’ 

Here, we examine the microscopic origin of the bi- 
collinear magnetic order in FeTe. In previous calcula¬ 
tions based on the DFT 6,7 ', although it was suggested 
that the origin of the AFB order is attributed to large 
farther-neighbor magnetic interactions such as the third- 
neighbor magnetic interaction J 3 , it is unclear what in¬ 
duces the large farther-neighbor magnetic interactions. 
One possible origin of such large farther-neighbor inter¬ 
actions is the superexchange interaction that is induced 
by farther-neighbor hoppings. To examine the effects of 
farther-neighbor hoppings, we artificially cut off farther- 
neighbor hoppings. As shown in Fig. 3(a), we find that 
the AFB becomes the ground state regardless of the 
range of such hoppings. The robustness of the AFB indi¬ 
cates that the superexchange interactions induced by the 
farther-neighbor hoppings play minor roles in stabilizing 
the AFB order. 

Another possible origin of the farther-neighbor mag¬ 
netic interactions is Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interactions, which are mainly induced by the 
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Fig. 3. (Color online) Ground-state energy per site of AF in the 
mVMC for FeTe with the cGW-SIC parameter (a) as a function of 
cutoff length of t and (b) as a function of ratio for J, A j (J —> A j-J) 
(in meV). The nearest-neighbor, next-nearest-neighbor, and 4th- 
neighbor hoppings are included in ”nn”, ”nnn”, and ”4th” calcula¬ 
tions, respectively. Energy is measured from that of the AFG. The 
system size is N s = 8 x 8. 


Hund’s rule coupling. To examine the roles of the Hund’s 
rule coupling as well as of the RKKY interactions, we 
artificially change the magnitude of the Hund’s rule cou¬ 
pling by introducing the scaling parameter A j, where A j 
scales on-site exchange interactions as J —> A jJ. Fig¬ 
ure 3(b) shows the ground-state energy per site of sev¬ 
eral AF orders in FeTe as a function of A j. The AFB 
becomes more stable as J becomes larger. This result in¬ 
dicates that the RKKY-type interactions play key roles 
in stabilizing the AFB order. Since scattering decreases 
with the magnetic ordering, the resistivity decreases well 
below the magnetic transition temperature 41 ) particu¬ 
larly along the (tt/2,—tt/2) ferromagnetic line. 42,43 ) The 
structural transition 44 ) would enhance the energy stabi¬ 
lization through the changes in transfer parameters. The 
couplings between itinerant electrons and localized spins 
resulting in the RKKY interactions were examined on 
phenomenological levels. 45-47 ) 

Now, we discuss the magnetic properties of FeSe. As 
in the case of FeTe, the AFS phase becomes the ground 
state in the LDA calculation for FeSe, as shown in Fig. 4. 
For FeSe, we only show the results of the largest sys¬ 
tem size ( N s = 8x8), because the size dependence of 
the energy is not smooth in FeSe, while we estimate the 
thermodynamic limit to be not far from the result for 
TVs = 8x8. By eliminating the double counting of the 


self-energy, we find that the AFH phase has the low¬ 
est energy, which is consistent with previous ab initio 
calculations. 48,49 ) What is more characteristic of FeSe 
is that the AFH and AFB phases win the AFS state 
and they have nearly the same energy, i.e., the difference 
in the energy between the AFH and AFB is on the or¬ 
der of 10 meV and the energies of the AFS and AFG 
phases are also close within 25 meV. Such degeneracies 
were not properly captured in the previous DFT-GGA 
studies, 48,49 ) where the energies of the AFB and AFG 
phases are much higher. Because of the limitation of the 
system sizes and the scattered size dependence of en¬ 
ergy for FeSe, the energy stability of FeSe has a large 
uncertainty with an energy scale of 10 meV. Thus, this 
degeneracy does not contradict the experimental results 
that the magnetic order is not observed despite its large 
local moment, since it is probable that the short-range 
order of the AFB and AFH (and AFS as well) coexist as 
a mixed or fluctuating phase, which does not allow the 
development of the long-range order above the supercon¬ 
ducting transition temperature. 

It is naively expected that the superposition of the 
degenerate magnetic orders will induce spin-glass-like 
states. However, the spin-glass-like states are not con¬ 
sistent with the existence of small but clear Fermi sur¬ 
faces 12,13 ) In the magnetically degenerate region, it is 
plausible that an unconventional superconducting phase 
becomes stable. By examining the stability of the su¬ 
perconducting phase by extending the present analysis 
for the normal states, the origin of small Fermi surfaces 
and the resultant BEC-like superconductivity 14 ) will be 
clarified. Further theoretical investigation of such a pos¬ 
sibility is left for future study. It is also interesting to 
determine the shapes of the Fermi surfaces from jumps 
in momentum distributions by performing the calcula¬ 
tions for larger system sizes. 

The energy of the AFG phase is relatively low in FeSe 
compared with that in FeTe. The reason why the energy 
of the AFG phase is low is that the occupation num¬ 
ber of the XY orbital having a very weak frustration is 
also close to half filling (1.10) in addition to those of the 
YZjZX and X 2 — Y 2 having stronger frustrations. Note 
that the AFG order is more stable if the unfrustrated or¬ 
bitals are closer to half filling. This reduction of the XY 
filling originates from the remarkable decrease in the on¬ 
site potential of the 3 Z 2 — R 2 orbital in FeSe; thus, the 
filling of the other orbitals must decrease. However, the 
YZjZX and X 2 — Y 2 orbitals are already pinned near 
half filling. This drives the filling of the XY orbital closer 
to half filling. Because the XY orbital partially resolves 
the frustration in comparison to FeTe and the energy of 
the AFG phase becomes lower, FeSe shows an accidental 
degeneracy of several phases. More precisely, when the 
three orbitals ( YZ , ZX , and X 2 — Y 2 ) could become 
close to half filling by an appropriate on-site potential, 
the AFS, AFB, and AFH phases may have similar ener¬ 
gies. If the XY orbital in addition could become close to 
half filling, the AFG phase joins in this degeneracy. 

As we will show later, the degeneracy may be lifted 
easily under pressure. The difference in the crystal field 
splitting from FeTe yields the unique energetic degener- 
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Fig. 4. (Color online) Ground-state energy per site of AF in 
the mVMC for FeSe with f LDA and with ^GW-SIC (in meV). The 
left-hand data is obtained by using the cGW-SIC calculation, but 
the on-site potential of the 3Z 2 — R 2 orbital is increased by 100 
meV artificially. The system size is N a = 8x8. 


acy of the AF orders in FeSe. 

Lastly, we examine the robustness of the degeneracy 
in FeSe. By comparing the results of the LDA and cGW- 
SIC models, it is speculated that the stability of the AF 
phases largely depends on the on-site potential of the 
3 Z 2 — R 2 orbital in iron chalcogenides. To determine 
the roles of the on-site potential directly, we monitor 
the effect of the on-site potential of the 3 Z 2 — R 2 in 
the cGW-SIC model for FeSe. As a result, we find that 
the degeneracy is easily lifted and that the AFS becomes 
the ground state as shown in Fig. 4. This result can be 
understood by considering the fact that the occupation 
number of the YZjZX orbitals deviates from the half fill¬ 
ing, and that only the X 2 — Y 2 orbital, which stabilizes 
the AFS states through large next-nearest transfer, stays 
near half filling and becomes the sole driving force of the 
magnetism. This mechanism is hardly captured by weak- 
coupling theories. 50 ^ Although not so easy, experimental 
control of the one-body potential of the 3 Z 2 — R 2 orbital 
will reveal the origin of the degeneracy in FeSe. The mag¬ 
netic transition under pressure 51,52 ) would be related to 
the degeneracy of the AF phases. In the present ab initio 
study using the cGW-SIC model of FeSe under 4 GPa, 
such a degeneracy is lifted and it predicts that the AFH 
phase becomes stable, while the AFS phase is competing 
with it. 53 ) 

In this Letter, we have shown that the elimination of 
the double counting in the exchange-correlation energy 
is essential for reproducing the AFB order observed in 
FeTe. The stability of the AFB in FeTe is insensitive 
to farther-neighbor transfer but sensitive to the strength 
of the Hund’s rule coupling. These results indicate that 
the nesting in Fermi surfaces and superexchange inter¬ 
actions play minor roles but that the RKKY interac¬ 
tions induced by the Hund’s rule coupling play a key 
role in stabilizing the AFB order. We have also found 
that several magnetic orders are energetically degener¬ 
ate in FeSe. This degeneracy is the possible origin of the 
absence of magnetic order in FeSe. It is shown that this 


degeneracy can be lifted by changing the chemical po¬ 
tential of the 3 Z 2 — R 2 orbital. Our analysis provides 
a firm theoretical basis for understanding the character¬ 
istic magnetism in iron chalcogenides and offers a clue 
to understanding other unique behaviors observed in the 
iron chalcogenides such as the possible high-temperature 
superconductivity in a single layer of FeSe. 15 ) 
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FeSe 



Fig. S 1. (Color online) Electronic band structures of the Fe 3 d 
MLWFs of FeSe in the LDA [(red) solid line] and the cGW-SIC 
[(green) dashed line]. The Fermi energy is set to zero. 



Fig. S 2. (Color online) On-site potential of Wannier orbitals for 
FeSe and FeTe obtained by the LDA and cGW-SIC calculations. 


local-density-approximation (LDA) calculation. We also 
show the transfer integral of the MLWFs in the con¬ 
strained GW with the self-interaction correction (cGW- 
SIC) model in Table S IV. Figure S 1 is the correspond¬ 
ing band structure of FeSe. The on-site potential of each 
model is shown in Fig. S 2. The detail of the cGW method 
is explained in Refs. SI and S4. 


FeSe 



Fig. S 3. (Color online) Ground state energy per site of AF in 
the mVMC for FeSe at ambient pressure and at 4.0 GPa with 
jcGW-SIC ( in me y) System size is N s = 8 X 8. 


S.2 Magnetism of FeSe under Pressure 

We calculate the magnetism of FeSe under pressure as 
well as that at ambient pressure. Figure S 3 shows the 
energy per site by the cGW-SIC calculation for FeSe at 
4.0 GPa in the experimental geometry. S5 ) Although the 
AFS-type short-ranged fluctuations were observed S6,S7 ) 
at ambient pressure, the unique energetic degeneracy at 
ambient pressure is lifted under pressure, and we predict 
that the AFH phase becomes slightly more stable than 
the AFS. The experimental verification would be desired. 


Supplemental Materials 

S.l Values of Parameters of Low-Energy Effec¬ 
tive Model 

In this supplemental material, we show the values 
of the parameters of the ab initio low-energy effective 
model. Table S I (II) is the three-dimensional effective in¬ 
teraction of the maximally localized Wannier functions 
(MLWFs) originating from the Fe 3 d orbitals of FeSe 
(FeTe). S1 -* We obtain the two-dimensional effective in¬ 
teraction by uniformly subtracting 0.6 (0.4) eV from the 
three-dimensional effective interaction of FeSe (FeTe) fol¬ 
lowing Refs. S2 and S3. Table S III is the transfer inte¬ 
gral of the MLWFs of FeSe and FeTe obtained by the 
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Table S I. Bare and effective three-dimensional Coulomb interactions between two electrons for all the combinations of Fe 3d orbitals 
in FeSe (in eV). Here, v and J v represent the bare on-site and exchange Coulomb interactions, respectively. The static limit of the 
effective on-site and exchange Coulomb interactions are denoted by U and J, respectively. 


FeSe 

V 

U 


XY 

YZ 

3Z 2 - R 2 

ZX 

X 2 - Y 2 

XY 

YZ 

3Z 2 - R 2 

ZX 

X 2 - Y 2 

AY 

18.65 

16.50 

17.28 

16.50 

16.51 

4.51 

3.19 

3.21 

3.19 

3.48 

YZ 

16.50 

16.97 

17.05 

15.70 

15.30 

3.19 

4.11 

3.53 

3.02 

2.98 

3 Z 2 - R 2 

17.28 

17.05 

19.09 

17.05 

16.00 

3.21 

3.53 

4.68 

3.53 

3.01 

ZX 

16.50 

15.70 

17.05 

16.97 

15.30 

3.19 

3.02 

3.53 

4.11 

2.98 

x 2 - y 2 

16.51 

15.30 

16.00 

15.30 

15.88 

3.48 

2.98 

3.01 

2.98 

3.77 


Jv 

J 


XY 

YZ 

3Z 2 - R 2 

ZX 

X 2 - Y 2 

XY 

YZ 

3Z 2 - R 2 

ZX 

X 2 - Y 2 

XY 


0.66 

0.79 

0.66 

0.34 


0.57 

0.69 

0.57 

0.32 

YZ 

0.66 


0.46 

0.56 

0.61 

0.57 


0.42 

0.48 

0.53 

3 Z 2 - R 2 

0.79 

0.46 


0.46 

0.75 

0.69 

0.42 


0.42 

0.62 

ZX 

0.58 

0.56 

0.46 


0.61 

0.57 

0.48 

0.42 


0.53 

X 2 - Y 2 

0.34 

0.61 

0.75 

0.61 


0.32 

0.53 

0.62 

0.53 



Table S II. Bare and effective three-dimensional Coulomb interactions between two electrons for all the combinations of Fe 3 d orbitals 
in FeTe (in eV). Here, v and J v represent the bare on-site and exchange Coulomb interactions, respectively. The static limit of the 
effective on-site and exchange Coulomb interactions are denoted by U and J, respectively. 


FeTe 

V 

U 


XY 

YZ 

3Z 2 - R 2 

ZX 

A' 2 - Y 2 

XY 

YZ 

3Z 2 - R 2 

ZA 

A 2 - Y 2 

XY 

17.08 

15.05 

16.20 

15.05 

16.26 

3.46 

2.30 

2.36 

2.30 

2.80 

YZ 

15.05 

15.36 

15.87 

14.24 

14.94 

2.30 

3.08 

2.62 

2.15 

2.30 

3Z 2 - R 2 

16.20 

15.87 

18.25 

15.87 

16.08 

2.36 

2.62 

3.73 

2.62 

2.35 

ZX 

15.05 

14.24 

15.87 

15.36 

14.94 

2.30 

2.15 

2.62 

3.08 

2.30 

X 2 - Y 2 

16.32 

14.94 

16.08 

14.94 

16.77 

2.82 

2.30 

2.36 

2.30 

3.39 


Jv 

J 


XY 

YZ 

3Z 2 - R 2 

ZX 

A 2 - Y 2 

XY 

YZ 

3Z 2 - R 2 

ZA 

A 2 - Y 2 

XY 


0.59 

0.73 

0.59 

0.33 


0.49 

0.62 

0.49 

0.31 

YZ 

0.59 


0.42 

0.49 

0.59 

0.49 


0.37 

0.40 

0.49 

3Z 2 - R 2 

0.73 

0.42 


0.42 

0.74 

0.62 

0.37 


0.37 

0.62 

ZX 

0.59 

0.49 

0.42 


0.59 

0.49 

0.40 

0.37 


0.49 

>■< 

1 

<N 

0.33 

0.59 

0.74 

0.59 


0.31 

0.49 

0.62 

0.49 
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Table S III. Transfer integral for the 3d orbitals of the Fe sites in the FeSe and FeTe, t]^^(Rx, Ry , Rz), where i LDA is the expectation 
value of the Kohn-Sham Hamiltonian for the Wannier function : t LDA = (</>|'H LDA |</>), m and n denote symmetries of the 3d orbitals, 
and the axis of ( R x ,R y ,R z ) is taken along the Fe-Se/Te directions. Units are given in meV. 


FeSe 


(m, ri) \ R 

[0,0,0] 

[§>4> o] 

[1,0,0] 

[1. 

-i,o] 

[I. 

T.°] 

[o, o, f] 

[1- 

-I £l 

2’ a J 

cry 

i 

a L 

(XY,XY) 

-509 

-410 

-70 


-11 


3 

-25 


6 

+ 

+ 

+ 

(XY, YZ) 

0 

273 

131 


-9 


-6 

0 


-9 

+ 

- 

-(1,4) 

(XY,3Z 2 -R 2 ) 

0 

-347 

0 


22 


-8 

0 


11 

- 

+ 

+ 

(XY, ZX) 

0 

273 

0 


-9 


18 

0 


-3 

- 

- 

“(1,2) 

(XY,X 2 -Y 2 ) 

0 

0 

0 


0 


-9 

0 


-4 

- 

+ 

- 

(YZ, YZ) 

46 

197 

128 


-17 


-8 

8 


27 

+ 

+ 

(4,4) 

(YZ,3Z 2 -R 2 ) 

0 

-119 

0 


7 


2 

0 


11 

- 

- 

-(4,3) 

(YZ, ZX) 

0 

127 

0 


-23 


-19 

0 


12 

- 

+ 

(4,2) 

(YZ,X 2 -Y 2 ) 

0 

223 

0 


1 


-3 

0 


20 

- 

- 

(4,5) 

(3Z 2 — R 2 ,3Z 2 — R 2 ) 

-388 

-4 

-15 


-14 


-6 

-23 


-9 

+ 

+ 

+ 

(3Z 2 -R 2 ,ZX) 

0 

119 

199 


-7 


-13 

0 


-10 

+ 

- 

“(3,2) 

(3 Z 2 ~ R 2 ,X 2 -Y 2 ) 

0 

0 

-115 


0 


1 

-8 


-6 

+ 

+ 

- 

(ZX, ZX) 

46 

197 

335 


-17 


13 

8 


0 

+ 

+ 

(2,2) 

(ZX,X 2 -Y 2 ) 

0 

-223 

82 


-1 


-15 

0 


7 

+ 

- 

(2,5) 

(X 2 -Y 2 ,X 2 -Y 2 ) 

-34 

-56 

93 


0 


17 

-28 


4 

+ 

+ 

+ 

FeTe 


(m, n) \ R 

[0,0,0] 


[1,0,0] 

[1. 

-i,o] 

[I. 

-i«] 

[0,0, |] 

[h 

-I £l 

2’ aJ 

g y 

i 

a L 

(XY,XY) 

-452 

-378 

-11 


-41 


-1 

-31 


12 

+ 

+ 

+ 

(XY, YZ) 

0 

237 

109 


3 


-6 

0 


-12 

+ 

- 

-(1,4) 

(X Y,3Z 2 -R 2 ) 

0 

-336 

0 


33 


-9 

0 


21 

- 

+ 

+ 

(XY, ZX) 

0 

237 

0 


3 


35 

0 


-1 

- 

- 

-(1,2) 

(XY,X 2 -Y 2 ) 

0 

0 

0 


0 


-14 

0 


3 

- 

+ 

- 

(YZ, YZ) 

44 

156 

103 


-15 


-12 

13 


37 

+ 

+ 

(4,4) 

(YZ,3Z 2 -R 2 ) 

0 

-122 

0 


17 


6 

0 


11 

- 

- 

-(4,3) 

(YZ, ZX) 

0 

101 

0 


-27 


-25 

0 


14 

- 

+ 

(4,2) 

(YZ,X 2 — Y 2 ) 

0 

178 

0 


0 


-10 

0 


22 

- 

- 

(4,5) 

(3Z 2 - R 2 ,3Z 2 - R 2 ) 

-480 

-73 

-53 


3 


8 

-67 


-23 

+ 

+ 

+ 

(3 Z 2 - R 2 , ZX) 

0 

122 

198 


-17 


-17 

0 


-32 

+ 

- 

-(3,2) 

(3 Z 2 - R 2 ,X 2 -Y 2 ) 

0 

0 

-29 


0 


-6 

30 


-30 

+ 

+ 

- 

(ZX, ZX) 

44 

156 

300 


-15 


42 

13 


9 

+ 

4- 

(2,2) 

(ZX,X 2 -Y 2 ) 

0 

-178 

136 


0 


-23 

0 


26 

+ 

- 

(2,5) 

(X 2 -Y 2 ,X 2 -Y 2 ) 

-211 

66 

52 


9 


12 

16 


-24 

+ 

+ 

+ 
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Table S IV. Transfer integral for the 3d orbitals of the Fe sites in the FeSe and FeTe, t < ^^ v ~ slc (R x , Ry, Rz), where £ cGW " SIC i s the 
transfer integral without double-counting : £ cG W-SIC _ ^|^ C GW-SIC|^ ? m an d n d eno te symmetries of the 3d orbitals, and the axis of 
( R x , Ry, R Z ) is taken along the Fe-Se/Te directions. Units are given in meV. 


FeSe 


(m, n) \ R 

[0,0,0] 

[l.-l.o] 

[1,0,0] 

[1. 

-i,o] 

[I. 

T.°] 

[o.o. S] 

[I. 

-I £l 

2’ a J 

cry 

i 

a L 

(XY, XY) 

-751 

-466 

-20 


24 


3 

-38 


6 

+ 

+ 

+ 

(XY, YZ) 

0 

201 

106 


-26 


-3 

0 


-9 

+ 

- 

-(1.4) 

(XY,3Z 2 -R 2 ) 

0 

-391 

0 


27 


-3 

0 


8 

- 

+ 

+ 

(XY, ZX) 

0 

201 

0 


-26 


-9 

0 


-4 

- 

- 

“(1,2) 

(XY,X 2 -Y 2 ) 

0 

0 

0 


0 


6 

0 


-10 

- 

+ 

- 

(YZ, YZ) 

-137 

97 

148 


-35 


-15 

7 


27 

+ 

+ 

(4,4) 

(YZ,3Z 2 -R 2 ) 

0 

-94 

0 


22 


-1 

0 


8 

- 

- 

-(4,3) 

(YZ, ZX) 

0 

199 

0 


-37 


-36 

0 


16 

- 

+ 

(4,2) 

(YZ, X 2 — Y 2 ) 

0 

204 

0 


-1 


7 

0 


17 

- 

- 

(4,5) 

(3Z 2 — R 2 ,3Z 2 — R 2 ) 

-1075 

-67 

-76 


5 


11 

-4 


3 

+ 

+ 

+ 

(3 Z 2 -R 2 ,ZX) 

0 

94 

195 


-22 


-16 

0 


-3 

+ 

- 

“(3,2) 

(3Z 2 - R 2 ,X 2 -Y 2 ) 

0 

0 

-35 


0 


-13 

-22 


7 

+ 

+ 

- 

(ZX, ZX) 

-137 

97 

263 


-35 


25 

7 


-2 

+ 

+ 

(2,2) 

(ZX,X 2 -Y 2 ) 

0 

-204 

147 


1 


-39 

0 


0 

+ 

- 

(2,5) 

(X 2 -Y 2 ,X 2 -Y 2 ) 

-299 

186 

19 


-20 


8 

-35 


17 

+ 

+ 

+ 

FeTe 


(m, n) \ R 

[0,0,0] 

[§,-§> °] 

[1,0,0] 

[1. 

-i,o] 

[I. 

T.°] 

[0.0, f] 

[i 

-I £l 

2’ a J 

Gy 

i 

a L 

(XY,XY) 

-678 

-410 

52 


5 


-2 

23 


-13 

+ 

+ 

+ 

(XY, YZ) 

0 

133 

79 


-11 


-1 

0 


11 

+ 

- 

-(1,4) 

(XY,3Z 2 -R 2 ) 

0 

-340 

0 


26 


-8 

0 


11 

- 

+ 

+ 

(XY, ZX) 

0 

133 

0 


-11 


-2 

-26 


10 

- 

- 

“(1,2) 

(XY,X 2 -Y 2 ) 

0 

0 

0 


0 


-1 

0 


16 

- 

+ 

- 

(YZ, YZ) 

-146 

42 

120 


-28 


-17 

11 


15 

+ 

+ 

(4,4) 

(YZ,3Z 2 -R 2 ) 

0 

-88 

0 


24 


6 

-23 


14 

- 

- 

-(4,3) 

(YZ, ZX) 

0 

186 

0 


-35 


-46 

0 


16 

- 

+ 

(4,2) 

(YZ,X 2 — Y 2 ) 

0 

139 

0 


1 


-2 

12 


14 

- 

- 

(4,5) 

(3Z 2 -R 2 ,3Z 2 -R 2 ) 

-974 

-114 

-114 


21 


26 

-1 


-1 

+ 

+ 

+ 

(3 Z 2 - R 2 , ZX) 

0 

88 

181 


-24 


-12 

0 


-4 

+ 

- 

-(3,2) 

(3 Z 2 - R 2 ,X 2 -Y 2 ) 

0 

0 

49 


0 


-24 

-1 


2 

+ 

+ 

- 

(ZX, ZX) 

-146 

42 

213 


-28 


51 

22 


13 

+ 

+ 

(2,2) 

(ZX,X 2 -Y 2 ) 

0 

-139 

147 


-1 


-30 

0 


8 

+ 

- 

(2,5) 

(X 2 -Y 2 ,X 2 -Y 2 ) 

-495 

292 

-47 


-9 


20 

5 


-13 

+ 

+ 

+ 
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